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ABSTRACT 

In  catalysis  theory  there  is  interest  in  the  number  of  solutions  of  the 
equation 

v"  +  A<1  +  $  -  v)Pe“Y/V  , 
with  the  boundary  conditions 

v' (0)  -  0,  v<1)  -  1  , 

the  parameters  A,  8,  Y  being  all  positive  and  p  a  non-negative  integer. 
The  paper  answers  this  question  when  Y  is  large,  which  is  the  interesting 
situation  physically,  and  although  the  treatment  is  somewhat  different  in  the 
cases  p  "  0  and  p  *  0,  the  final  answer  is  the  same,  that,  given  8, 
there  exist  two  positive  functions  A^  (y)  and  A^(y)  such  that  the  problem 
has  one  solution  if  A  <  A^(y)  or  A  >  three  solutions  if 

A,j(Y)  <  A  <  A2(y),  and  two  solutions  if  A  *  A^y)  or  A  *  A^ty). 

AMS  (MOS)  Subject  Classifications:  34E99,  80A20,  80A30 
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SIGNIFICANCE  AND  EXPIJtNATION 

In  the  production  of  chemical* ,  catalysts  are  often  required  to  convert 
gaseous  reactants  into  useful  products.  Frequently  the  catalyst  is  in  the 
font  of  a  porous  pellet  and  the  gas  must  diffuse  into  the  interior  of  the 
pellet  so  that  the  catalyst  there  is  fully  utilized.  Depending  upon  the 
relative  rates  of  diffusion  and  reaction,  temperature  and  concentration 
gradients  are  set  up  across  the  pellet,  and  their  determination  is  essential 
for  the  calculation  of  the  over-all  rate  of  conversion,  The  modelling  of 
these  processes  within  the  pellet  leads  to  a  set  of  parabolic  partial 
differential  equations,  and  a  first  step  in  the  study  of  these  is  to  determine 
whether  there  exist  steady-state  solutions,  and,  if  so,  how  many  of  these 
there  are. 

The  present  paper  works  at  a  particular  one-dimensional  steady-state 
equation  which  nonetheless  seems  to  be  typical  of  more  general  situations,  and 
it  is  shown  rigorously  that  if  the  activation  energy  is  sufficiently  high, 
then  the  number  of  solutions  must  be  essentially  either  one  or  three 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sussury  lies  with  HRC,  and  not  with  the  authors  of  this  report. 
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THE  NUMBER  OF  SOLUTIONS  OF  AN  EQUATION  FROM  CATALYSIS 


S.  P.  Hastings  and  J.  B.  McLaod 


1.  INTRODUCTION 

Na  ara  intarastad  in  tha  nusbar  of  solutions  of  the  aquation 
(1.1)  v*  ♦  X{1  +  8  -  v)Pa"Y/T  , 


with  tha  boundary  conditions 

(1.2)  v'(0)  -  0,  v(1)  -  1  . 


This  is  tha  one-dimensional  casa  of  an  aquation  of  soon  importance  in  catalysis  theory ,  tha 
parasMtars  1.  8,  Y  baing  all  positive,  and  p  a  non-negative  integer.  Readers 
intarastad  in  tha  derivation  of  the  aquation  and  its  physical  interpretation  ara  referred 
to  the  work  by  Aria  (1),  particularly  section  2.S.4. 

Relevant  work  on  tha  analysis  of  this  problem  has  bean  dona  by  Dancer  (2)  and  Barter 
(3].  Dancer’s  work  is  concerned  with  tha  casa  p  •  1,  although  ha  mentions  that  it  can  be 
extended  to  other  values  of  p,  and  he  shows  that,  with  8  fixed  and  y  large,  there  are 
at  moat  three  solutions  of  the  problem  except  for  a  relatively  ssuill  range  of  values  of  X 
where  he  is  unable  to  coaw  to  any  conclusion.  His  argument,  which  depends  upon  regarding 
(1.1)  as  a  perturbation  of  the  Gelfand  equation 

u*  +  UeU  -  0 

and  using  ideas  from  bifurcation  theory,  extends  also  to  radially  symmetric  solutions  in 
two  dlmnsions,  and  he  shows  that  the  problem  is  fundamentally  different  in  higher 
dimensions  still,  in  that,  for  y  large,  there  may  be  values  of  X  giving  large  numbers 
of  solutions.  His  work  in  one  and  two  dimensions  is  therefore  a  partial  answer  to  the 
conjecture  that  in  these  dimensions  the  problem  has  at  most  three  solutions. 

Barter  is  concerned  with  two  dimensions  and  an  equation  which  is  a  generalisation  of 
the  case  p  «  0.  His  results  delineate  a  region  of  the  parameter  space  in  which  there 
exist  at  least  three  solutions,  and  another  in  which  there  exists  at  most  one. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  and  in  part  by 
Grant  No.  DAJA37-81-C-0220.  Also  supported  by  the  U.  K.  Science  and  Bigineering  Research 
Council. 


Our  results  are  more  siailar  to  those  of  Dancer,  although  obtained  in  quite  a 
different  way.  The  restriction  to  one  diaension  arises  because  we  intend  to  use  the 
autonoaous  nature  of  the  equation  in  that  case  to  reduce  the  problea  to  an  integration,  and 
we  can  then  establish  that  for  (1,1)  with  p  ■  0  and  the  boundary  data 
(1.3)  v* (0)  -  0,  v(1>  -  0  , 

the  problea  has  either  two  solutions  or  none,  this  regardless  of  whether  or  not  Y  is 
large.  (Given  Y,  there  is  one  exceptional  value  of  1  for  which  the  two  solutions 
becoae  coincident.)  The  introduction  of  the  boundary  condition  v(1)  -  1  can  be  regarded, 
at  least  for  large  Y,  as  a  perturbation,  and  by  studying  this  perturbation  we  are  able  to 
assert  that  (1.1)-(1.2)  has,  if  Y  is  sufficiently  large,  either  one  solution  or  three 
(except  for  two  exceptional  values  of  X  where  there  are  two).  Unlike  the  Dancer  result, 
there  is  no  restriction  on  X.  There  would  seen  little  doubt  that  the  result  raulna  true 
whether  Y  is  large  or  not,  but  the  analysis  required  for  its  proof  seems  too  formidable. 

Our  approach  is  similar  in  spirit  to  that  of  Smoller  and  Wasserman  [4],  who  studied 
the  number  of  solutions  of  another  autonomous  problem  by  reducing  it  to  an  integral.  Their 
nonlinearity  was  cubic,  however,  resulting  in  completely  different  analysis. 

For  p  *  0,  and  boundary  data  (1.3),  we  establish  again  that  there  are  either  two 
solutions  or  none,  although,  in  contrast  with  the  case  p  -  0,  we  have  to  assisse  Y  large 
in  order  to  assert  this.  This  in  turn  leads  to  the  same  result  as  for  p  -  0  with  the 
original  boundary  conditions  (1.2),  namely  that  for  Y  large  there  is  either  one  solution 
or  three,  except  at  two  transitional  values  of  X  where  there  are  two. 

There  is  a  variant  of  the  problem,  still  of  physical  interest,  in  which  the  equation 
is  replaced  by 

v*  ♦  X { 0  +  ^  *  v(1)(l  -  -  v}pe  -  0  , 

with  the  boundary  conditions 

v*(0)  -  0,  v'(1)  -  P { 1  -  v( 1 )}  , 

U,  v  positive.  This  has  been  studied  by  Kapila  and  Matkowsky  (5],  and  also  by  D.  G. 
Schaeffer  (private  communication] ,  who  show  that  the  perturbation  corresponding  to  changing 
from  (1.3)  to  (1.2)  can  produce  two  further  solutions,  so  that  there  can  be  from  one  to 
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five.  Thia  problaai  ia  not  atudied  in  tha  praaant  pa  par,  although  thara  would  aeon  to  ba  no 
raaaon  why  tha  aathoda  uaad  hara  ahould  not  axtand . 

-V  /y 

tha  arguaenta  would  eartainly  allow  ua  to  altar  tha  nonlinaarlty  a  to  aoaathing 

-v/v 

■ora  general,  hut  thara  doaa  aaaa  to  ba  a  larga  aaaaura  of  agreement  that  a  ia  tha 
correct  nonlinaarlty,  and  ao  no  attempt  ha  a  baan  made  to  push  tha  thaoraaa  in  tha  direction 
of  ganarallaation. 

In  $2  wo  prove  tha  baaic  theoraa  for  (1.1)  and  (1.3)  with  p  -  0,  and  carry  out  tha 
perturbation  to  (1.2)  in  {3.  Later  aectione  daal  with  tha  caaa  p  >  0. 
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2.  TH»  CAM  p  ■  0  WITH  ZERO  BOUNDARY  DATA 

*i  are  interested  in  the  boundary-value  problem 

•y/v 

(2.1)  V  *■  Xa  -  0  , 

(2.2)  v'(0)  -  0,  v(1)  -  0  . 

Sine*  X  >  0,  and  tha  exponential  ia  non-negative,  we  aniat  have  v*  <  0,  and  ao  v'  <  0, 
v(0 )  >  0.  Sat 

v  «  aY  , 

where  a  ia  ao  choaen  that 

Y(0)  -  1  . 


The  equation  for  Y  ia 


Set 


and  we  have 


niua 


a  H  ♦  >«'T/“V  -  0 

d*2 


1/2 

x  -  a  X,  Y*  -  dY/<3X  , 


y-  .  x«'Y/aY  -  0  , 

Y(0)  -  1,  y*(0)  -  0,  Y <a“1/2>  -  0 


r 


-2X  /  e’Y,/atdt  , 


f  --  <2Xf1/2(J  .*^dt)-,/2  . 

Y 


and  integration  with  respect  to  Y  over  [0,1]  gives 


(2.3) 


-1/2  -  (2X,-1/2  l  (/’  e-Y/otdt)-1/2dY  , 
0  Y 


or,  with  Y  “  «P, 
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<2. 4) 


(f)1/2  -  p-1/2  l'  U  .•g/tat)-,/2dr  . 

'  o  y 

(tow  our  object  ia  to  ahow  that,  given  X,  y ,  there  exiat  at  moat  two  aolutiona  a  of 
(2.3),  or,  equivalently,  at  Boat  two  aolutiona  U  of  (2.4),  for,  through  v(0)  •  a,  two 
aolutiona  for  a  a eana  two  aolutiona  for  our  original  boundary-value  problem.  What  we 
ehall  actually  prove,  and  it  clearly  iapliea  the  above,  ia 
Theorem  1.  For  |i  >  0,  define  f(M)  b^ 

f(u)  -  u- 1/2  f'  (/’  e-^dt )-1/2dY  . 

0  Y 

Then  f*(u:  haa  preclaely  one  aero,  and  that  elmple. 

Proof.  Let  P(u)  “  e“u  .  Then 


(2.5) 


df 

U  m 
*1 


"I*  2 


/  J-P’^dt 

±«'/2l  Jh - 

0  1/  P(£)*t}3/2 


dY 


Since  F'  «  -F,  and  t  <  1  in  the  rangea  of  integration,  we  see  that 


(2.6) 


J  F(£)dt «  -  /  J  p’(*r)<*t  , 


and  so 


Hence  f'(U)  >0  if  u  >  1, 


«  >  .  1  t  +  1  f  . 
dp  2  2 


and  we  need  concern  ouraelvea  only  with  u  <  1. 


Then 
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<*  i* 


i  4L  I  „V2  f  JL 

-  2  “  «  "  I  M  I  1 


1  [ 


dY 


{/  *($«*! 

Y 


3/2 


,  I/'  ^F-(^dt}2  ,  /  ^ 

3  .  3/2  f  _Y _  -  i  u3/2  f  -£~ -  <W 


+  i  „3/2  /  -  dY  -  j  v'  I 

4  0  If1  .(EWlV2  0  If  k(E)  dt} 3/2 


U  *(*)<*}' 

Y 


But 


(2.7) 


(2.8) 


Y  t 


/  1_  r.(E)at  -  -  i  lr-(»)  -  <  -  l  {»Mb)  -  "•($)}  , 

/’  £  r*(£)dt  -  -  i  (nu)  -  **(*)}  ♦  l  J  r(*)at  , 


so  that 


and 


T  -4  »"({9«t  <  -  {'  i  r*(£)at  *  -  /  , 


Y  t 


,  J 


dY 


,  I/'  i  K“)«l2  , 

3  y 3/2  |  Jf - _  dy  ♦  1  „3/2  J  ”4 - - 

4  •  I)’  0  1/ 


dY  - 


When 


—  «  0,  we  can  use  (2.5)  to  obtain 


dy 


Vt 
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i  /  £'•(*)* 

-  i  uV2  /  -Y-  —  ■  -  dv 

0  1/  *(iw\3/2 

y 


i  «/  ?*’£)*)*  ,  /  i  »*(»)* 

1UV2/  - <»,  .3  nj'.l\  J  dY 

°  1/  <)<*}5/2  0  {/  P(£)at}3/2 

Y  ¥ 


>  (}  »'/2  +  I  U3/2  -  »3/2 )!  -V - — 

0  {/  r(*)dt}3/2 


tram  (2.6).  Hence,  with  w  <  1,  f'(P)  -  0  lmpliaa  ftp)  >  0,  and  so  f'(ti)  has  at 

■ost  ona  zero,  and  that  simple. 

It  remains  to  be  proved  that  f  does  possess  one  zero.  However  for  small  p  we 


f(D)  -  2u"1/2,  f(P>  ~  -m‘3/2  , 

and  we  have  already  seen  that  f'(p)  >  0  for  p  >  1.  so  that  certainly  f'(p)  has  a 


The  following  result  Is  an  Immediate  corollary  of  Theorem  1. 

Theorem  2.  For  the  boundary- value  problem  (2.1)- (2.2),  there  exists  a  positive  number  K 

X/y  <  K,  Just  one  solution  if  XA  “  K.  end  two 


solutions  if  X/Y  >  K. 


3.  THE  CASE 


0  WITH  NON-ZERO  BOUNDARY  DATA 


Ha  now  consider 


v*  +  Xe“Y/w  -  0  , 


v'(0)  -  0,  v(1)  -  1  . 


Hith  tha  sue  substitutions  as  in  {2,  this  reduces  to 

Y-  ♦  Xe-Y/0Y  -  0  . 


Y(0)  -  1,  Y'(0)  -  0,  Y(a-1/2)  -  a”1,  a  >  1  , 


which  in  turn  leads  to 


or,  with  Y  •  ay,  to 


Let  us  now  define 


<2lf1/2  J1  [/'  e-Y/atdt)-1/2dY  , 


(^),/2  -  u-1/2  j'  {/’  .-u/tdt}-1/2dY  . 

u/y  y 


f,(U)  -U-1/2  /’  U  e-U/tdtrV2dY  . 

UA  Y 

He  are  interested  only  in  Y  large,  and  then  f  (u>  behaves  like  ftp)  in  Theorem  1 
except  when  p  also  is  large.  In  particular,  there  1b  an  e  >  0  such  that,  if  Y  is 
sufficiently  large,  then  on  the  Interval  0  <  p  <  CY,  f  1  (p)  first  decreases  to  a  minimum, 
and  then  increases.  On  the  other  hand,  while  f ( p)  is  defined  for  all  p  >  0,  with 
f(p)  ♦  •  as  p  ♦  ■>,  f^p)  makes  sense  only  for  U  <  Y»  and  f^Y)  “  0.  Our  goal  is  to 
show  that  (3. 1)-(3.2)  has  at  most  three  solutions,  as  stated  in  Theorem  4  below,  and  this 
is  a  clear  consequence  of 

Theorem  3.  For  Y  sufficiently  large.  f^(p)  has  precisely  two  taros  in  the  range 
0  <  P  <  Y,  and  those  are  simple.  One  occurs  close  to  the  zero  of  f'(u),  and  the  other 


P  ”  Y  -  2  log  Xq  +  0 (Y  log  Y>  , 

where  x0  is  the  unique  root  exceeding  1  o£ 
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(3.3) 


cosh”'x  -  x/(x2  -  1)1/2  . 


(That  there  ie  a  unique  root  exceeding  1  of  (3.3)  followe  by  considering  the  difference 
of  the  two  aides/  recalling  that 

f-  cosh_1x  -  (x2  -  1)“1/2  .) 
ax 


Theorem  4. 

If 

Y  is  sufficiently  large  in  the 

boundary-value  problem  (3.1)— (3.2),  then 

there  exist 

two  positive  functions  K.(Y),  JC„(Y) 

such  that  the  problem  has  one  solution  if 

X  <  K,(Y) 

or 

-  l  A 

X  >  K. (Y),  three  solutions  if 

K. < Y )  <  x  <  (Y),  and  two  solutions  if 

X  -  K^Y) 

or 

X  -  Kj(Y). 

Proof  of  Theorem  3. 

Making  the  changes  of  variable 

t  •  u*.  Y  “  UZ,  c  -  u"\  Z  -  o-1  , 

we  have 


(3.4) 


1/U  VM 

f  {/ 

1/Y  Z 


T  t  (u )  «  J  {/  e  ^'dt}  ,/,2dZ  , 


and 


f'/v  -t/x  „-v., 

J  e  d*  «  e  J  e  ds 


-U  r  1  U-U  .  -u  r 

e  J—  e  du  •  e  J 


u-v  -6 


d6 


U  u 


0  (9  +  U) 


Now,  by  integration  by  parts. 


u-u 


-0 


d0 


(0  +  U) 


,  U-u 

.  S - } 

V  U2  J 


U-u 


-0 


2  J 
0 


(0  +  u> 


d0 


so  that 


(3.5) 


U-M  -9  U-u 

J  - - - r  d9 { l  +  r(u,U)>  -  . 

0  ( 0  +  U )  u  u 


where 
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f(p,U)  -2  / 


U-p  -0 


0  (0  +  p ) ' 


d0  /  / 


U-p  -6 


0  <6  +  p ) 


It  is  clear  that 


f(p,U)  <  2p  , 


but  we  can  also  assert  that,  for  large  u, 

(3.7)  (p.U)  -  0(U~2)  . 

uniformly  in  U  for  p  <  U  <  Y.  (He  recall  that  we  are  only  interested  in  large  p, 
since  f^(p)  behaves  like  ftp)  when  p  is  not  large.)  The  truth  of  (3.7)  is  obvious  so 
far  as  those  terms  in  the  differentiation  are  concerned  which  arose  from  differentiation 
under  the  integral  signs.  The  terms  arising  from  differentiating  the  limits  of  Integration 
are  (omitting  a  factor  2) 


p-U  U-p  -0  p-U  U-p  -0  U-p  -0 

1-V'  — t-T^+s-rJ  — — J  -«»/{/  — £— : 

U  0  (0  +  p)  U  0  (0  +  p)  0  (0  +  p) 


p-U  U-P  -0  U-p  -0 

V  J  (0  -  p  -  0>d0  /  {/  —S—j 

U  0  (0  +  p)  o  (0  +  p) 


If  U  -  p  <  1 ,  then  the  numerator  does  not  exceed 


p-6  /  (U-p  -  0)d8  -  —  p  (U-p)  , 

0 

while  the  denominator  is  not  less  then 

p-U  U-p  ,  .  , 

(S-Z~  /  d0 ) 2  >  Xp  4(U  -  p)2 
U  0 

for  some  constant  K,  since  U  -  p  <  1.  Thus  (3.7)  follows  if  U  -  p  <  1.  If 


U  -  p  >  1 ,  the  numerator  does  not  exceed 


-6  P-U  r  -8,.  _  -1  -6 

p  (U  -  p)eH  /  e  d6  <  e  p  , 

0 


while  the  denominator  is  not  less  than 


1  -8 


u  —* l— 7 


0  (8  +  p) 


(P  ♦  1> 
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Hence  again  (3.7)  follows 


from  (3.4),  (3.5),  (3.6), 


TV  Y  ,  ,  U-U  . 

f  <u)  -  e  /  ~  (  -  *-r-j  /2{i  ♦  0(v  >>d0 

U  0  V  0 


1  _1  “U 

ly  ~rf-Tv75  {1  +  0<“'1>)i 

P  (0«  -ps)' 


”  ,  2--24p-i/l^>0(p-^,)ax, 

1  (X  -  p  s  ) 


whsrs  x  “  Ue*  and  ths  O-term,  thanks  to  (3.7),  has  the  property  that.  If  differentiated 

1 

-2  F 

with  respect  to  p,  it  yields  o(p  ).  Thus,  with  x  -  pe  t. 


f,(P)  «  e2  / 

'  e 


1  7  F 

F  ,Y®  /w®  2uo“2  -1 

s2  /  .....m.  -  {,  ♦  0(p  ’))dt 

t  <t  -  i)  ' 


If  p  »  y  -  log  y«  then  0  -  p  *  0{log  p),  and  so 


f  (P)  -  2p~  e  / 

*  4 


x  7 .  7 

.  Ye  /pe  ,  ... 

“V  /  (t^  -  1)  Wi{x  +  0(p  log  p ) }dt 


-1  F  -1  F  F  -i 

2p  e  cosh  (Ye  /pe  ){1  +  0(p  log  p)}  , 


our  previous  renarka  about  differentiating  the  O-term  implying  now  that  differentiation  of 


the  last  O-term  will  yield  0(p  log  p).  Thus 


♦  OCU-1  )}f ,(u> 


.  2_y 


1 


♦  2u  1,2  1/2(-  ^fr^®2  ^2  M  +  ’){1  +  o(i’"1 109  w,)  • 


w  e 


u  e 


and  ao  f'(u)  -  0  whan 


-i  r  71 

coah  (Ye  /lie  ) 


2  , 


2  Y  1/2 

(^-  1) 


{1  +  0<y  1  log  y )}  , 


2  y 

U  a 


l.e.  whan 


2 1  ?  -1 
Ye  /lie  •  xQ  +  Oty  log  y)  , 

whare  xQ  la  daflned  In  tha  atatamant  of  the  theorem.  Thus 

U  -  Y  -  2  log  Xq  ♦  0(Y-1  log  Y)  • 

Me  are  now  assured  that  there  la  a  root  of  f *  C u )  with  y  >  y  “  log  Y»  and  that  it 
necessarily  satisfies  (3.9).  To  prove  that  there  is  just  one,  and  that  one  simple,  we 
evaluate  f^(y)  with  y  satisfying  (3.9)  and  show  that  f*(y)  <  0.  (We  leave  the 
straightforward  calculations  to  the  reader.)  He  can  also  use  (3.8)  in  similar  manner  to 
prove  that  f  ^ ( U )  >0  if  ey  <  u  <  y  -  log  Y»  and  the  proof  of  Theorem  3  is  complete. 
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4.  TH1  CA81  p  *  0  WITH  ZIRO  BOUNDARY  DATA 


Ms  ar*  Interested  in  ths  problaa 

v«  +  1(1  ♦  0  -  v)pe‘Y/v  , 


with 


vMOJ  -  0,  v(1)  -  0  , 

and  v#  ar*  lntsrsstsd  in  solutions  (or  which  v  <  1  +  0.  As  in  $2,  ws  sst 

v  -  ay,  x  -  a1/Sx  , 


whsrs 


T<0)  -  1,  a  -  v(0)  <  1  . 

Ths  equation  for  Y  (Y*  “  dY/dX)  is 

Y*  ♦  i(t  ♦  #  ■  aY)ps"Y/aY  -  0  , 

Y(0)  -  1,  YMO)  -  0,  Y(a'V2>  -  0  . 

Thus 

a"1/2  -  (2l)"1/2  /  (/  (1  ♦  B  -  at)ps"Y/atdt)‘1/2dY  , 

0  Y 

or,  with  r  -  aw,  a  <  1  ♦  B,  w  >  Y/<1  ♦  B), 


(4.2) 


-  W 


-1/2 


1  1 

/  (/ 

0  Y 


(1  +  B  -  Yt/W)Ps-|,/tdt]1/2dY  . 


Given  X ,y ,  ws  want  to  show  that  thsrs  exist  at  most  two  solutions  a  of  (4.1),  or  at 
most  two  solutions  w  of  (4.2).  Mo  shall  actually  prove 
Theorem  5.  For  w  >  Y/d  *  B)  define  g(w)  b£ 

g(U)  -  w"V2  /  (/  (1  +  B  -  Yt/w)Ps”M/tdt)“1/2dY  . 

0  Y 

Then  for  Y  sufficiently  large,  B  fixed,  g' (w)  has  precisely  one  sero.  and  that 
simple.  Further,  for  this  sero,  u  is  close  to  Y/( 1  *  B),  being  given  mors  precisely  by 
(4.9)  or  (4,10)  below. 

Proof.  For  convenience,  write 

K  -  (1  ♦  B)A,  t  -  W/u,  Y  -  w/T  , 


1J- 


and  we  have 


With 


I 

P 

Y  g(w) 


/’  f2{JT  u'2(k 

u  u 


i)P.-“du}-1/2 


dr  . 


u  »  0  ♦  v,  T  •  a  *  v  , 


m  have 


Yg(U)  *  /  <U  ♦  o)’2{/  (U  ♦  8)”2(k  -  ■--?  g)Pa~|J~fld8}~1'/2<lg 
0  0  wo 


1  1 

r  •  / 

o 


(w  ♦  of2{/ 


HL 


8  -  K~ 1 ) P 


a  At 


P+2 


.-V172 


do 


(4.3) 


G(W)  - 


F 


(w  +  o) 


-2U 


UL  »  .8  -  K 


V  -e 


(W  +  8) 


P+2 


ae} 


-1/2 


do 


1 

F 

where  G(u)  “  (KY)  g(M). 

Wa  shall  show  that  G’ (w)  (or  equivalently  g'(w))  can  ba  aero  only  if  vi  -  K  1  ia 
aaall.  and  wa  first  Investigate  the  inner  integral  in  (4.3).  For  convenience,  we  write 


and  than  the  inner  integral  ia 


(4.4) 


/ 


0 


(W  ♦  8)P*2 


-8 

e 


d8 


»*  / 
a 


a+o 


p  -T 

t*e 


(u  +  t  -  a) 


P+2 


dT  . 


To  deal  first  with  a  large,  we  note  that  if  we  differentiate  G(m),  we  can 
1 

F  -2 

differentiate  e  (which  Merely  Multiplies  by  1/2),  or  we  can  differentiate  (w  ♦  o) 
(which  Multiplies  by  a  factor  0(u  'll,  or  we  can  differentiate  the  inner  integral. 
Differentiating  the  left-hand  expression  in  (4.4)  under  the  integral  sign,  we  see  that  we 
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f 


collect  a  factor  0(a  1 )  or  0(p  ^),  end  so  whan  we  differentiate  G(p)  we  similarly 

collect  euch  a  factor,  and  eo  the  leading  tare  in  the  derivative  la  that  arlaing  from  the 
1 

y 

differentiation  of  e  ,  and  G* (p)  >  0. 

To  deal  next  with  the  caae  where  a  la  moderate,  we  see  that 

(  (4.4)  St  0app‘<pf2>  If  a  <  A  , 


(4.S) 


(4.4)e“*  X  !  ^e^dt  If  o  >  A  . 

a 


where  A  la  a  fixed  number  chosen  so  that 

(a  ♦  o)pe  *  for  c  >  A  , 

and  X  swans  that  each  aide  la  bounded  by  a  positive  constant  tlaiea  the  other  as  u  ♦ 

Me  remark  also  that  If  we  differentiate  (4.4)  (with  raapect  to  p ) ,  then  we  see  that  the 
order  relatlona  obtained  by  tonally  differentiating  (4.S)  are  valid.  (This  uses  the 
choice  of  A.)  Thus 


G(w)  *  eh^V  *  '  (M  ♦  c,-VV2dc  ♦  K-'U  rP.-TdT)-1/2} 
0  a 


*  y\*u  tp.-'«.r,/2 . 


In  view  of  the  remark  just  made  about  differentiation,  we  can  differentiate  this  last  order 
relation,  whence  it  is  clear  that  G' (p )  >  0. 

Me  know  therefore  that  a  must  be  small  for  G' (p)  •  0.  To  detenine  a  more 
precisely,  we  write 


G(p) 


tJ 

0 


4  I  ' 


-  G1(p)  +  G  (p)  ♦  Gj(M)  , 


say,  where  e  is  a  small  fixed  positive  number,  and  A  a  large  one,  both  independent  of 
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ii.  Than 


ip-1  c 

G  ^  (u )  »  a2  y2  {t-KXu-1  >}{1+0(e))  ]  j  -  —  -  [{•♦o)p+1-  apf1l|"1^i4» 

0  P 


(p*1)1/2a2'lw^P'  a2  **  f  f  ((♦+1)p*1-  l}"1/2a#J1>0lw“1)}{l-KMc)} 
0 


The  integral  in  the  laat  formula  line  convargea,  to  Cp,  nay,  aa  It  ♦  •  and  a  ♦  0,  e 
being  fixed,  provided  that  p  >  1.  If  p  “  1,  it  ia  asymptotic  to 

log  t  -  log  a  . 

It  ia  easy  to  verify  that  formal  differentiation  of  the  aaymptotic  expreeaion  for 
Gf  ia  juetified,  and  that  tha  major  contribution  to  Gj  cornea  from  the  differentiation  of 
factor a  involving  a.  Thua,  for  p  >  1, 


1  1/2  2”  F*1  -  2<P+tl  £/* 

G*  <|i)  »  t  (p  ♦  1)  '  (1  “  p)a  v  a 
1  *  0 


I  {(♦  +  D1*1  -  x 


"1/2  , 


*  {1  +  Oty*’)}(l  ♦  0(£>){1  ♦  0(a)) 


1/2  F  F-1  -  2(p41’  2 

-  (p  ♦  1 >v  V  F  *  (e/a)  *  (e/a  )  x 


x  {1  +  0(it-1)){1  +  o(e>){1  ♦  o( a/e ) )  , 


and  again,  aa  u  ♦  •  and  a  ♦  0,  e  fixed,  the  major  term  in  G*(u)  above  ia  the 


firat.  for  p  »  1, 


1  1  1 

1/2  F  F"  ’  -1  1/2  F  -1/2  -1 

G’(u)  -  -(p  +  1)ve  u  a  -  -2 v  a*  u  '  a 

for  g2,  there  exiat  conatanta  X^(e,A),  Kj(e,A)  auch  that 

11.  11. 

F  F"1  F  F-1 

e  v  x, (e ,h)  <  G2(u)  <  a  u  K2(e,A)  , 


11.  11. 

F  F"1  F  F~ 

a  It  K^e.A)  <  Gj(ll)  <  •  U  KjU.A) 


-16- 


Clearly ,  G'^  daminatea  Gj. 

For  G,,  o  >  A,  and  ao 


(4.6) 


(4.7) 


v 

A 
1-1 


G,(n)  “a2  /  (w  ♦  o)  2{/  — e”9d»}  1^2do 
3  *  o  (w  ♦  e )**2 


-  .a  f  <m  *  or’uf  -  n  At----i2  a-^#»d.}',/2d 

A  0  0  (|l  ♦  8)*^ 


k-'  - 


iv* 


-  •  I  (v  *  «>■*{(/  ~  ]  ) 

A  a  a+o  <u  -  a  +  t)*^ 


V1  - 


A  a+o  a  ^  -t 


■*  /  (M  ♦  Of2!!/  ♦  /  -  /  )  - E*  _ 

A  0  A  0  (|l  -  a  ♦  T)*^ 


dt}"1/2dO 

dTj-’^dO 


1^,-1  1  t 

a2  m2*^  /  (U>«)"2(/  TPa'Tdt)'1/2da(140{u‘,)){la0(a)Kia0(APa“*)}  , 

a  0 


where  the  0-terma  aay  depend  on  A, 

1-1  1 


•(pi 


— p 

)",/2e2  ii  (1  ♦  0(m*')){1  ♦  0(a)){1  ♦  0(APe"A))  . 


In  differentiating  GJ(  we  must  be  careful,  for  formal  differentiation  of  the  leading  tare 
doee  not  give  the  correct  reeult.  If  we  denote  the  expreeelon  in  {....)  in  (4.6)  by 
<1(0, u),  then 


V1  - 


*(M)  ■  -  2e2  /  (ii  ♦  o)-3(J(o.M)}-,/2do 


(4.6) 


1-1 


-  i  e2  f  (M  ♦  o)~2(j(«.i,)}-i/2|',-i-y>P«^W>  -  ££.}  do 

A  (u  «•  o,^2  II1’*2 


In  paaalng  froai  (4.6)  to  (4.7),  we  aew  that,  for  large  it.  A, 

,-(p*2> 


J(o,u)  ~  pin 


1-1  1 


2  jr 

and  ao  the  firat  term  on  the  right  of  (4.8)  la  0(e  v  ),  which  cannot  balance 
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G’ (u) •  The  term  involving  (a  +  o)P  is  again 


1-1  3 


1-1  1 


o(e2  u2  /  (u  o)”^p*^ope”°do)  -  o(a2  u2**"  )  , 


but  tha  tarn  involving  ap  ia  asymptotically 

{  £  ,(pd"3/2»’p*V  (» . -  j  .  ?)-' . 


For  GMu)  >0  we  must  hava  a  balanca  between  G|(|t)  and  G£(u)»  and  so,  fixing  a,  A 
as  small  as  w  please,  and  than  latting  u  ♦  a  ♦  0,  we  hava,  for  p  >  1, 

J  (P  +  1)V2(p  -  Da2'* u2**'  a  2<Pf  )  I  {(♦  +  1)pf1  -  1}"1/2<H 

0 

if1  *- 


-1 


1  ,  ..-3/2  2*  P  D 
~  r  (pl )  a  U  r  , 


so  that 


(4.9)  a2<  **  ’  ~  (pl)3/2(p  -  IHp  ♦  l)1^-1-—^  /  ((♦  *■  D5*1  -  1>-1/2d*  . 

T  0 


tn  tha  caaa  p  •  1,  wa  hava 


1  1  , 

1/3  ?,  ?"  _-1 


1-1  1 


(p  ♦  1)  a  w  a 


1  ,  ,,-3/2  2*  P  p 
~  -  (pi)  a  u  a*^  , 


a  ~  23/4«1  ♦  #)VVV2 


ao  that 
(4.10) 

Finally,  to  show  that  tha  caro  of  G’(u)  is  uniqua  and  simpla,  wa  hava  to  vsrify 
that,  in  the  relevant  range  for  u,  G*(l>)  >  0.  This  is  now  a  routine  calculation,  and  wa 
leave  it  to  tha  reader. 
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G*<U> 


P 


Y-U 


<U  +  o) 


-2 


{/ 

0 


(U  ♦  9  -  K 


-1 


<U  ♦  9) 


P+2 


^  .-ed9|-1/2do 


Theorem  6.  For  Y  sufficiently  largo,  6  fixed,  K  «  (1  +  B)/y»  G*'(u)  hae  preclaely 

two  zeros  In  the  range  Y/M  ♦  8)  <  u  <  Y,  both  simple.  One  occurs  near  the  zero  of 
G(u)»  the  other  asymptotically  at  U  “  Y  -  2  log  xQ,  where  Xg  le  as  In  Theorem  3. 
Proof.  The  proof  of  this  need  not  be  given  in  detail  since  it  is  in  effect  just  a 
repetition  of  parts  of  the  proofs  of  nieoreas  3  and  5,  If  u  is  such  that  Y  -  I*  ♦  •  as 
Y  ♦  •»  then  G*(u)  can  be  analysed  as  G(u)  in  Theorem  S.  The  result  is  thet,  in  such 
a  range  of  u,  G*'(u)  has  precisely  one  zero,  and  that  staple,  with  the  value  of 
a  •  u  -  K  1  satisfying  <4. 9)  or  (4.10).  If  u  is  such  that  Y  -  U  “  o(Y)»  and  this 
clearly  overlaps  with  the  previous  range  of  u,  then  we  can  treat  G*  as  we  did  f  ^  in 
Theorem  3.  What  the  analysis  there  shows  is  that,  to  determine  the  asymptotic  behaviour  to 
the  required  degree  of  accuracy,  we  treat  u  ♦  9  -  K  '  as  though  it  were  y  -  X  \  and 
U  +  9  and  y  +  o  as  though  they  were  u.  Then 


1  1 


GNU)  ~  FF'iv  -  K-V  ?  r  U  «"8d0)"V2dO 


1  1  , 
P  P-1 

e  u 


0 


(u  -  X-1)  1  I  (1  -  e-0)-1/2do 


1  1  « 

P  p-1  -1  '  P  , 

2e  m  <y  -  X  )  2  / 

1 


1  exP4<Y-U>) 


(,2  -  „-1/2dx 


1  1  1 

P  p-1  -1  ‘  P  -1  P  P 

-  2«  u  (u  -  X  )  cosh  (e  /e  )  . 


This  is  the  iaM  as  the  asyaptotic  expression  for  f.(y)  except  that  u  is  replaced  by 
1  1  11  11 
P*1  -i  "  P  V  P  V  P 

u  (y  -  X  )  ,  and  Ye  /ue  by  e  /e  The  latter  pair  are  asymptotically  the 
same  when  Y-U  is  bounded,  which  is  the  range  we  are  interested  in,  and  the  former  pair 
play  no  part  in  determining  the  asymptotic  position  of  u  since  when  we  differentiate 
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these  factor*  In  f|  or  G* ,  thair  darivativaa  do  not  contribute  to  the 

leading  terms  in  f  ^  or  G* ' .  Thus  the  final  answer  for  G*  is  asymptotically  the  same 

as  that  for  fj,  as  required. 
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